Evidence for exhaustion in the conductivity of the infinite-dimensional periodic 

Anderson model 



A.N. Tahvildar-Zadeh\ M. Jarrell\ Th. Pruschke^ and J. K. Freericks^ 
^ Department of Physics, University of Cincinnati, Cincinnati, Ohio 
^ Institut fiir Theoretische Physik, Universitdt Regensburg, Regensburg, Germany 
Department of Physics, Georgetown University, Washington, DC 20057 
(February 1, 2008) 



Monte Carlo-maximum entropy calculations of the con- 
ductivity of the infinite-dimensional periodic Anderson model 
are presented. We show that the optical conductivity dis- 
plays anomalies associated with the exhaustion of conduction- 
band states near the Fermi energy including a Drude weight 
which rises with temperature, with weight transferred from a 
temperature and doping-dependent mid-IR peak and a low- 
frequency incoherent contribution. Both the Drude and mid- 
IR peaks persist up to very high temperatures. The resistiv- 
ity displays a non-universal peak and two other regions as- 
sociated with impurity-like physics at high temperatures and 
Fermi Liquid formation at low T. 



I. INTRODUCTION 

Metallic compounds containing rare earth elements 
with partially filled / shells, such as CeBeia or UPta, be- 
long to the general category of heavy fermions materials 

. They are characterized by a large Pauli susceptibility 
and specific heat as compared to ordinary metals, which 
indicate a large effective electronic mass, and also by 
anomalous transport properties such as non-monotonic 
temperature dependence of the resistivity. These anoma- 
lies are usually attributed to the formation of a reso- 
nant state at the Fermi energy, associated with moment 
screening, due to the admixture of the electronically ac- 
tive and well localized / orbitals with the metallic band of 
the host. As the Kondo screening peak develops, the op- 
tical conductivity develops a pronounced low-frequency 
Drude peak. Interband transitions across the Fermi en- 
ergy can also yield strongly temperature-dependent fea- 
tures in the mid-IR optical response. 

Heavy fermion (HF) materials are usually modeled by 
the periodic Anderson model (PAM) or the single impu- 
rity Anderson model (SIAM). While the SIAM is able 
to capture the general physics of moment formation and 
subsequent Kondo screening and is thus sufficient to un- 
derstand most of the thermodynamics of HF materials, 
lattice effects like coherence and correlations between dif- 
ferent sites are of course completely neglected in such a 
simplified picture. Coherence effects lead for example 
to a vanishing low-temperature resistivity in the metal- 
lic regime and to the formation of a band insulator in 



the symmetric PAM. Correlations between the sites on 
the other hand are responsible for various types of phase 
transitions which lead to an extremely rich phase dia- 
gram with metallic, insulating or near insulating regimes 
depending upon the filling and model parameters. 

In the metallic regime of the PAM, perhaps the most 
important but also controversial lattice effect is the so 
called exhaustion as proposed by Nozieres. The "exhaus- 
tion" problem |^,^ occurs when a few mobile electrons. 
User, have to screen many local moments, n^, in a metal- 
lic environment, i.e, Uscr << nf. This situation is en- 
gendered by the fact that only the electrons within Tk 
(where Tr- is the single impurity Kondo temperature) of 
the Fermi surface can effectively participate in screening 
the local moments. Thus the number of screening elec- 
trons can be estimated as Uscr — Nd{Q)TK, where Nd{0) 
is the conduction band density of states at the Fermi 
level. A measure of exhaustion is the dimensionless ratio 
p — — jv^(c/)Tk Hi' Nozieres has argued that p is 
roughly the number of scattering events between a local 
moment and a mobile electron necessary for the mobile 
electron's spin to precess around a local moment by 27r. 

In the case p S> 1, when the number of screening elec- 
trons is much smaller than the number of local moments 
to screen, magnetic screening is necessarily collective and 
the single impurity picture becomes invalid. 

One further consequence of Nozieres' idea then is Q| 
that the small fraction riser of screened states may be 
viewed as relatively stable polaron-like particles at low 
temperatures T < Tk- He thus proposed that the 
screened and unscreened sites may be mapped onto parti- 
cles and holes of an effective single-band Hubbard model, 
respectively. The screening clouds can hop from site to 
site and may effectively screen all the moments in a dy- 
namical fashion. 

In previous work, we extended [Q this argument to 
explain the strong reduction of the Kondo scale ||] in 
the metallic regime of the PAM. The hopping constant 
of Nozieres' effective model is suppressed relative to the 
bare one by the overlap of the screened and unscreened 
states. This, together with the fact that the filling of 
the effective Hubbard model is very close to half filling 
(rie// ~ rif — User ~ 1 for n/ = 1), suggests that the rel- 
evant low-energy scale Tq of the effective model becomes 
very small, compared to the Kondo impurity temperature 
Tk. 
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In the crossover regime characterized by temperatures 
between Tk, where conventional Kondo screening be- 
gins, and the lattice scale Tq, where coherence forms, 
the quasiparticle peak and the screened local moment 
evolve much more slowly than their counterparts in the 
SIAM [Q. This provides a possible explanation for the 
slow evolution of the quasiparticle peak seen in photoe- 
mission experiments However, these results remain 
controversial and are complicated by the fact that 
photoemission is very surface sensitive 1^, and probes 
only the first ~ 10 A of the surface under investigation. 
Thus, predictions for less surface sensitive probes such as 
transport and optical conductivity are important. 

In this paper we study the effect of screening on the 
optical conductivity and resistivity of the PAM. We 
find that the resistivity is characterized by an impurity- 
like regime at high temperatures T « Tk, a crossover 
regime with non-universal behavior, and a low tempera- 
ture regime T ^ Tq where coherence begins to form. The 
crossover regime is the most interesting and is charac- 
terized by the exhaustion of screening states near the 
Fermi surface, as measured by the development of a 
dip centered at the Fermi energy, w = 0, in the effec- 
tive f-d hybridization r(cL'). The peak in the resistiv- 
ity, which marks the upper end of the crossover regime, 
is a non- universal feature centered at a temperature Te, 
To < Tg < Tk- Te is the temperature at which exhaus- 
tion first becomes apparent as a reduction in r(0). In the 
exhaustion region, the Drude weight increases strongly 
with temperature, and the lost weight at low T is largely 
transfered to a mid infra-red (MIR) peak. This feature 
is due to interband transitions and its position matches 
with the minimum of the direct energy gap between the 
two quasi-particle bands below and above the Fermi en- 
ergy. Both the Drude and MIR peaks persist up to un- 
usually high temperatures T w lOTo. We interpret the 
persistence of these features as well as the unusual de- 
cline of the Drude weight as T — as manifestations of 
Nozieires' exhaustion scenario. 



neighbors and scaled as t = t*/2y'D. We take t* = 1 as 
our unit of energy. U is the screened on-site Coulomb re- 
pulsion for the localized / states and V is the hybridiza- 
tion between d and / states. This model retains the 
features of the impurity problem, including moment for- 
mation and screening, but is further complicated by the 
lattice effects. 

Metzner and VoUhardt observed that the irre- 
ducible self-energy and vertex-functions become purely 
local as the coordination number of the lattice increases. 
As a consequence, the solution of an interacting lattice 
model in D — oo may be mapped onto the solution of 
a local correlated impurity coupled to a self-consistently 
determined host We employ the quantum Monte 

Carlo (QMC) algorithm of Hirsch and Fye ||l2| to solve 
the remaining impurity problem and calculate the imagi- 
nary time local Green's functions. We then use the max- 
imum entropy method (MEM) to find the / and d 
density of states and the self-energy |14|. In the limit 
of Z? = oo the two-particle vertex corrections to the op- 
tical conductivity (ct^)) vanish identically due to sym- 
metry restrictions |lq|. Hence (t(w) can be calculated 
by the knowledge of one-particle self energy through a 
particle-hole bubble diagram which includes the appro- 
priate electron-photon interaction vertices evk at each 
end. A straightforward calculation of this diagram leads 

to ra, 



f{oj)-,f{e + uj) 
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(2) 



where V is the lattice volume, / is the Fermi function, 
Ad{ek,oj) — —-Im[Gd{k,uj)] is the conduction band 
spectral function and p{y) — exp{~y'^) / ^/^T is the non- 
interacting density of states. 



III. RESULTS 



II. FORMALISM 

The Hamiltonian of the PAM on a D-dimensional hy- 
percubic lattice is 

H = — ^ V ( d]d,„ + H.c.) 

i(7 1(7 

+ Ec^(«M-i/2)(n/a-i/2) . (1) 

i 

In Eq. (1), d{f)\^J destroys (creates) a d{f) electron with 
spin a on site i. The hopping is restricted to the nearest 



We calculated the transport for a variety of model pa- 
rameters. The QMC-MEM results are limited to rela- 
tively small values of U/V by the QMC procedure and 
relatively low temperatures where MEM may be used to 
calculate the dynamics. Here we present QMC results 
for U — 1.5, V = 0.6, « 1 for three conduction- 
band fillings Ud — 0.4, 0.6, and 0.8. The impurity scale 
Tk was calculated by extrapolating the local susceptibil- 
ity to T = (i.e. Tk = l/x™p(r 0)). The lattice 
Kondo temperature was similarly calculated by extrapo- 
lating the effective local susceptibility of the PAM which 
is the additional local susceptibility due to the introduc- 
tion of the effective impurity into a host of d-electrons 
0. These extrapolations result Tk = 0.2,0.46,0.5 and 
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To = 0.017,0.055,0.149 for the three band fiUings men- 
tioned above, respectively. 



A. Optical Conductivity 

The real part of the optical conductivity calculated 
with QMC-MEM is plotted versus frequency in Fig. |l| 
for several temperatures. As the temperature is lowered 
T < Tfc two features begin to develop, a Drude peak at 
zero frequency, associated with free quasiparticles, and 
a Mid-IR peak aX uj k, \ generally associated with an 
interband transition as we will discuss later. 




FIG. 1. Optical conductivity for various temperatures when 
U = 1.5, V = 0.6, n/ ^ I, Ud = 0.6 and To = 0.055. 
Both Drude and Mid-IR peaks are present and persist up 
T /To ^ 10. The inset shows the corresponding band energies 
Ek plotted vs. the bare d-band dispersion ek when T /To = 0.46. 
The solid line shows the quasiparticle energy calculated as the 
real part of the Green's functions poles. The symbols show the 
positions of the maxima in the f and d spectral functions. The 
mid-IR peak in the optical conductivity is due to direct tran- 
sitions between occupied and unoccupied quasiparticle states. 

To be consistent with the way that experimental data 
is usually analyzed the features in cr(u;) are fit to a 
Lorentzian plus two (asymmetric) harmonic-oscillator 
forms for the higher-energy peak 
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with T the relaxation time of the "free" quasiparticles, 
and the constants Cmir, ^mir, and Tmir are respec- 
tively the weight, center, and width of the mid-IR peak. 
As shown in Fig. ^(a), in order to improve the fitting 
procedure we include a second harmonic-oscillator form 
centered at ujinc ~ 0.1. Since this weight cannot be fit 
to the Drude form, part of the low frequency spectral 



weight remains incoherent as T is lowered. An alterna- 
tive method to separate the contributions of the Drude 
and the mid-IR peak to the spectral weight is to fit the 
optical conductivity data to a "generalized" Drude form 
and a harmonic-oscillator with strength C^jjj^ for the 
mid-IR region as discussed above 
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For a generic Fermi-liquid model a = 2 as was assumed in 
the previous fitting form. But, it turns out that treating 
a as a fitting parameter renders the use of incoherent 
harmonic-oscillator unnecessary as shown in Fig. ^. The 
inset to this figure shows the value of exponent a as a 
function of T. 
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FIG. 2. Two fitting schemes for separating the share of 
Drude and interband transition to the optical conductivity. 
The model parameters are the same as in the previous fig- 
ure and T /To — 0.46. a) Fitting to Eq. a Fermi-liquid 
Drude form and two harmonic-oscillator forms. One of the 
harmonic-oscillators accounts for the interband transitions in 
the mid-IR region and the other for the low frequency spec- 
tral weight that remains incoherent, b) Fitting to Eq. vA, a 
"generalized" Drude form for the low frequency part and a 
harmonic-oscillator for the mid-IR region. Using the expo- 
nent in the "generalized" Drude form as a fitting parameter 
(shown in the inset) renders the use of an incoherent har- 
monic-oscillator peak unnecessary. 

The spectral weights D^, Cmir, Cine, as well as those 
from the "generalized" fit D* and Cl^jj^ are plotted in 
Fig. H when U — 1.5, V = 0.6, n/ ~ 1 and = 0.6 
(To = 0.055). The most surprising result is that there 
is a significant transfer of weight from high to low fre- 
quencies with increasing temperature. The weight in the 
Drude peak, as measured by either Dq or D* , rises with 
temperature at low T. This feature was independent of 
the fitting procedure applied and was even found if the 
extrapolation technique of Scalapino et al. was em- 
ployed at low temperatures T < Tr-. The opposite trend 
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is generally expected since the Drude peak is associated 
with the buildup of quasiparticle weight in the single- 
particle spectra at low temperatures. 
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FIG. 3. The Drude, mid-IR and incoherent weights vs. 
temperature when U = 1.5, V = 0.6, Uf ~ 1, n^j « 0.6 
and To — 0.055. The coejficients Do, Cmir and dnc 
( D* and C'lufj) were determined from a fit to Eq, |^ 
(Eq. The inset shows the effective f-d hybridization 

T[lu) — Im (l/G/(ti^) -I- E(tj)) plotted vs. u> for several dif- 
ferent temperatures. The weight in the Drude peak falls with 
r(0) since the weight in the quasiparticle peak is expected to 
vary as ~ exp(-7rC//8r(0))/r(0) (see text). 

To interpret the drop in the Drude weight with low- 
ering T, we calculate the effective hybridization strength 



T{uj) = Im T,{uj) + 



1 



(5) 



where Gf{uj) is the local / Green's function and 
is the local f-electron self-energy. F(w) is a measure of 
the hybridization between the effective impurity in the 
DMF problem and its medium (for example, in the SIAM 
T{uj) = ttV'^ Nci{uj) , where Nci{llj) is the d-band density of 
states). The inset to Fig. || shows the effective hybridiza- 
tion for the PAM near the Fermi energy. r(aj) begins to 
develop a dip at the Fermi energy at roughly the same 
temperature where Drude weight begins to drop. 

The formation of this dip in T(lu) can be interpreted as 
an effective reduction or "exhaustion" of the states near 
the Fermi energy responsible for screening the local mo- 
ments. In infinite dimensions, other mechanisms which 
could be responsible for the dip, such as non-local spin 
or charge correlations, are absent |]l8|. We therefore be- 
lieve that it can be taken as direct evidence for Nozieres' 
exhaustion scenario. Concomitant with the development 
of this dip is a drop in the quasiparticle weight as may be 
seen from the following simple qualitative argument: The 
relevant low-energy scale is expected to vary roughly like 
the Kondo temperature of the effective impurity model 



with r(a;) as medium, i.e. Tq ~ exp (— 7rC//8r(0)). Since 
this scale also sets the width of the quasiparticle peak, 
whose height is on the other hand fixed to ^f^Xoy 
Friedel's sum rule, the weight in the quasi-particle peak 
is proportional to exp (— 7rC//8r(0)) /r(0), which falls ex- 
ponentially with decreasing r(0). Since the Drude peak 
is formed from quasiparticle excitations, the fall in the 
r(0) also corresponds to a loss in Drude weight. 

The weight lost in the Drude peak as T decreases is 
gained by both the mid-IR peak and the low-frequency 
incoherent part since we find the total spectral weight to 
be roughly constant in T. The weight gained in the inco- 
herent part is reflected by an increase in Cine and by the 
reduction of the exponent a. We observe that for T < Tp 
the exponent starts to deviate from its Fermi-liquid value 
2, indicating that additional weight is building up in the 
tail of the low frequency peak. The temperature depen- 
dence of D* is basically the same as that of Dq -\- Cine 
(cf. Fig. H) suggesting that there is a low-frequency in- 
coherent contribution to the spectral spectral function, 
in addition to the quasiparticle peak, due to non-Fermi- 
liquid nature of the problem in this parameter regime. 
The remaining transferred weight is gained by the peak 
in the mid-IR region which is due to excitations between 
the occupied and unoccupied bands states. Since the 
optical conductivity is constructed from a convolution 
of lattice spectral functions, the origin of this peak can 
be confirmed by inspecting the dispersion of these bands 
which is plotted in the inset to Fig. |l|. Here, the solid line 
is determined from the real-part of the poles of the lat- 
tice propagators, thus it represents the quasiparticle dis- 
persion. The location of the peaks in the single-particle 
spectra are denoted by the symbols. The singly occupied 
/-level appears as an almost dispersionless band below 
the Fermi energy. The energy of the mid-IR peak coin- 
cides with that of the interband transition between the 
two quasiparticle bands below and above the Fermi en- 
ergy at Efc ~ (this value of ek is special since it rep- 
resents most of the points in the Brillouin zone of the 
infinite dimensional lattice.) 

Both the Drude peak and the mid-IR peak in cr(w) 
persist up to unusually high temperatures T « IOTq. In 
the conventional Kondo picture, based upon the single- 
impurity model, the quasiparticle peak in the DOS, and 
hence both of these features in a{uj), are expected to dis- 
appear at much lower temperatures relative to the Kondo 
scale. The slow evolution of these features may be due 
to a crossover between two energy scales l^],^; Tr-, the 
screening scale of the impurity model with the same pa- 
rameters as the PAM, which is the temperature where 
screening starts, and Tq the Wilson-Kondo scale of the 
lattice. This screening regime is extended, over that in 
the SIAM, because of the reduction in r(0). I.e. as T is 
lowered below Tk , the Kondo scale of effective impurity 
problem in the DMFA is self-consistently reduced. Thus, 
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the slow evolution of Dq, Cmir and the quasiparticle 
peak itself are direct consequences of exhaustion. 

B. Resistivity 

The resistivity p = 1/ lim^_,o fT(w) also displays 
anomalies due to exhaustion, p versus T from a QMC- 
MEM calculation is shown in Fig. ^ in a semilog plot for 
three different conduction band fillings n^. 
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FIG. 4. The resistivity vs. temperature calculated with 
QMC-MEM when U = 1.5, V = 0.6, n/ ^ 1 for na = 0.4, 0.6 
and 0.8 (with To = 0.017, 0.055 and 0.149, respectively). The 
solid lines are from fits to a polynomial. At high temper- 
atures, the resistivity has an impurity-like log-linear regime 
which terminates maximum at T = Te. Te is non-universal 
in that it is independent of Tk or To which increase with Ud ; 
whereas Te decreases with Ud. When Ud = 0.4 and 0.6, Te 
corresponds to the temperature at which exhaustion becomes 
apparent as a dip in T{ll>). ForT < T^, the resistivity shows a 
non-universal crossover region (see text). When Ud = 0.4 and 
0.6, at still lower temperatures the resistivity begins to drop 
quickly atT ~ To, indicating the onset of Fermi-liquid forma- 
tion and coherence. However for Ud = 0.8 there is no such 
drop consistent with the proximity to the insulating regime 
(nj =nd = I). 

We can see the interplay between the different energy 
scales of the PAM in Fig. ^ and distinguish three distinct 
characteristic temperature regions for the resistivity, i) A 
region of log-linear resistivity which occurs at high tem- 
peratures T ^ Tk where the single impurity physics is 
dominant and the correlated sites act like independent 
(incoherent) scattering centers. As a result, the resistiv- 
ity displays the log-linear behavior of a dilute magnetic 
alloy, ii) As the temperature is lowered, the resistivity 
reaches a maximum value at T = Tg. Note that Tg is 
non-universal in that it is independent of Tk or Tq which 
increase with n^; whereas Tf, decreases with n^. At least 
for the Ud = 0.4 and 0.6 datasets, Tg corresponds to the 
temperature at which exhaustion first becomes apparent 
as a dip in As the temperature is lowered further, 



T < Te, a pronounced dip begins to develop in r(u;), and 
is accompanied by a rapid drop in the resistivity. Then, 
as the decline in r(0) slows, the resistivity falls more 
slowly. Since the impurity scattering rate is expected to 

vary like ~ In (^-^exp{-8U/nT{0))^ the re- 

sistivity can fall with decreasing r(0) even though the 
system has not begun to form a Fermi liguid. Partial 
Fermi-liquid formation can also contribute to the drop; 
however, this effect would be suppressed for small by 
exhaustion and the concurrent reduction of the effective 
Kondo scale. For Ud ~ 0.6 it is uncertain which of the 
two mechanisms has the dominant contribution to the 
drop in the resistivity at the top of the crossover regime. 
Also it is unclear why the crossover regime becomes less 
pronounced as Ud decreases. However, it seems likely 
that the drop for small n^, where exhaustion is most 
pronounced, is due the drop in the the scattering rate. 
Whereas for n^j — s- 1, where there is no exhaustion (the 
Kondo scale Tq is actually enhanced relative to the impu- 
rity scale Tk it seems likely that the initial drop 
is due to partial Fermi-liquid formation, iii) Finally, as 
the temperature is reduced further, for conduction band 
fillings Ud — 0.4 and 0.6, the resistivity begins to drop 
quickly at T ~ Tq, indicating the onset of Fermi- liquid 
formation at Tq for these fillings. However for Ud — 0.8 
there is no indication of this for temperatures well below 
Tq. When n/ = = 1, the PAM forms an insulating 
gap of width « Tq |0 (which is resonantly enhanced 
over the impurity scale Tk). Thus, for n^; < 1, we expect 
a drop in the resistivity to occur when the temperature 
equals roughly the energy difference between the chemi- 
cal potential and the bottom of the gap. 

IV. COMPARISON WITH EXPERIMENT 

Exhaustion should be most prevalent in heavy ferniion 
materials with low carrier concentration or a low density 
of conduction band states at the Fermi surface. One such 
material, for which an extensive study of the optical con- 
ductivity has been done, is Yb4As3 with a Kondo scale of 
To ~ 40K |2^. This material has strongly temperature- 
dependent Drude and MIR (at roughly 0.4 eV) peaks for 
temperatures ranging from 39K to 320K, as predicted by 
our model calculation. In addition, there is a very sig- 
nificant weight transfer from the Drude peak to the MIR 
peak as the temperature increases over this range 
There is also some evidence for a Drude weight which 
drops with lowering the temperature below T^ in a typi- 
cal paramagnetic HF material like CcAla Here, the 
transfer of weight is far less apparent, but may be in- 
ferred from the analysis of the data, where the inverse of 
the Drude weight (measured by fitting the optical con- 
ductivity to a Drude form at low frequencies) shows an 
increase with lowering T at w = 0. However, in this ma- 
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tcrial, the Drudc weight vanishes quickly as the temper- 
ature is increased to T = lOK « STq. Thus, exhaustion, 
and the concoinitant protracted evolution of D{T) would 
seem to be less evident in this material than it is in the 
low carrier system Yb4As3. Unfortunately, there is not 
enough optical data available for paramagnetic HF sys- 
tems to determine whether the drop in the Drude weight 
and other feature associated with exhaustion are generic 
features or even to determine their dependence on the 
f-levcl degeneracy and the conduction band filling. We 
believe that a systematic study of the optical properties 
of low carrier density HF systems could shed significant 
light on the importance of exhaustion in these materials. 

V. CONCLUSION 

We have performed a QMC-MEM simulation of the 

infinite-dimensional Periodic Anderson Model to study 
the optical conductivity in the metallic regime. We see 
several anomalies, including: (i) a Drude peak which 
persist up to anomalously high temperatures T < IOTq, 
(ii) a MIR peak due to interband transitions across the 
Fermi energy, which also persists up to high temperatures 
T < lOTo, (iii) a significant loss of Drude weight and f-d 
hybridization as T ^ 0. We interpreted these results in 
terms of Nozieres' exhaustion picture. Finally, we note 
that since optical measurements are far less surface sen- 
sitive than photoemission measurements, the presence of 
these features in experimental spectra could serve as less 
controversial evidence for exhaustion in heavy Fermion 
systems. 
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